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We calculate the two-loop pressure of a plasma of quarks and gluons at finite temperature and 
chemical potential using the hard thermal loop perturbation theory (HTLpt) reorganization of 
finite temperature/density quantum chromodynamics. The computation utilizes a high temperature 
expansion through fourth order in the ratio of the chemical potential to temperature. This allows 
us to reliably access the region of high temperature and small chemical potential. We compare 
our final result for the two-loop HTLpt pressure at finite temperature and chemical potential with 
perturbative quantum chromodynamics (QCD) calculations and available lattice QCD results. 
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I. INTRODUCTION 



Quantum chromodynamics (QCD) exhibits a rich phase structure and the equation of state (EOS) which describes 
the matter can be characterised by different degrees of freedom depending upon the temperature and the chemical 
potential. Hadrons are the relevant degrees of freedom at low temperature and chemical potential where chiral 
symmetry is spontaneously broken but the matter is SU(3)^ center-symmetric. At high temperatures the system 
is expected to make a phase transition to a quasifree state known as quark-gluon plasma (QGP). In the QGP 
chiral symmetry is restored and the center symmetry of SU(3)^ is spontaneously broken. At high temperatures 
and moderate chemical potentials one therefore expects the system to be in the QGP phase. Such conditions are 
generated in relativistic heavy ion collisions at Brookhaven National Laboratory's Relativistic Heavy Ion Collider 
(RHIC) [H, the European Organization for Nuclear Research's Large Hadron Collider (LHC) [H, and the Gesellschaft 
fur Schwerionenforschung's Facility for Antiproton and Ion Research (FAIR) [3|. 

The determination of the equation of state (EOS) of QCD matter is extremely important to QGP phenomenology. 
There are various effective models (see e.g. 04111) to describe the EOS of strongly interacting matter; however, one 
would prefer to utilize systematic first-principles QCD methods. The currently most reliable method for determining 
the EOS is lattice QCD At this point in time lattice calculations can be performed at arbitrary temperature, 
however, they are restricted to relatively small chemical potentials 0,0. Alternatively, perturbative QCD (pQCD) 
can be applied at high temperature and/or chemical potentials where the strong coupling {g^ = Airas) is small 
in magnitude and non-perturbative effects are expected to be small. However, due to infrared singularities in the 
gauge sector, the perturbative expansion of the finite-temperature and density QCD partition function breaks down 
at order requiring non-perturbative input albeit through a single numerically computable number [isl . [l6| . Up to 
order g^\w{l/g) it possible to calculate the necessary coefficients using analytic (resummed) perturbation theory. 

Since the advent ofpQCD there has been a tremendous effort to compute the pressure order by order in the weak 
coupling expansion [ij, E, 113, [13 ■ The pressure has been calculated to order of g^\D.{l/g) at zero chemical potential 
(/i = 0) and finite temperature T [l5| and finite chemical potential/temperature (/x > and T > 0) [l7| . In addition, 
the pressure is known to order g'^ for large ^ and arbitrary T (Tsj . Unfortunately, one finds that as successive perturba- 
tive orders are included, the series converges poorly and the dependence on the renormalization scale increases rather 
than decreases. The resulting perturbative series only becomes convergent at very high temperature (T ^ 10'^ Tc). 
One could be tempted to say that this is due to the largeness of the QCD coupling constant at realistic temperatures; 
however, in practice one finds that the relevant small quantity is, in fact, ag/Tr which for phcnomcnologically relevant 
temperatures is on the order of one-tenth. Instead, one finds that the coefficients of ag/Ti are large. This can be seen 
by examining the weak coupling expansion of the free energy J^{T, //) of QGP calculated [13 up to order ln(Q;s) 
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FIG. 1. The Nf = 3 pQCD pressure specified in Eq. IH) as a function of the temperature. Successive perturbative approxima- 
tions are shown through order agin as for vanishing fi (left) and for non- vanishing fi (right). The shaded bands indicate the 
variation of the pressure as the MS renormalisation scale is varied around a central value of A = I-k^T'^ + /i^/yr^ [13, [3 by a 
factor of two. We use A-j^ = 290 MeV based on recent lattice calculations [l^l of the three- loop running of as . 
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where here and throughout all hatted quantities are scaled by 27rT, e.g. /t ~ ii/{2-kT), A is the modified minimum 
subtraction (MS) renormalisation scale, and as = as (A) is the running coupling. At finite T the central value 
of the renormalisation scale is usually chosen to be 2'kT. Hovifever, at finite T and /i we use the central scale 
A = 27r-\/r2 _(_ (^/7r)2^ which is the geometric mean between 27rT and 2/x [13, [3- In Fig. [T] we plot the ratio of the 
pressure to an ideal gas of quarks and gluons. The figure clearly demonstrates the poor convergence of the naive 
perturbative series and the increasing sensitivity of the result to the renormalisation scale as successive orders in the 
weak coupling expansion are included. 

The slow convergence found when naively applying pQCD has led to various resummation schemes which attempt 
to improve the convergence of the successive approximations by reorganizing the calculation in terms of quasiparticle 
degrees of freedom [2ll - l37j . These resummation methods include some relevant physical ingredients, e.g. screening 
masses and Landau damping. These reorganisations of perturbation theory canonically include quasiparticle degrees 
of freedom from the outset, as opposed to naive perturbation theory. In the naive perturbative treatment an expansion 
around the vacuum is made and one only includes quasiparticle effects in order to regulate infrared divergenences. 
Based on Hard Thermal Loop (HTL) resummation (22l . [23j . a manifestly gauge- invariant reorganization of finite 
temperature/density QCD called HTL perturbation theory (HTLpt) has been developed [1^. HTLpt has so far been 
applied primarily to the case of finite temperature and zero chemical potential. In HTLpt [29l | the next-to-leading 
order (NLO) [lo] and next-to-next-to-leading order (NNLO) [3l| thermodynamic functions have been calculated at 
finite T but ^ = 0. Recently [H, the leading order (LO) HTL pressure for finite T and ^ has been calculated 
and approximately a decade ago it was applied at LO for finite /i but T = (37j . 

In view of the ongoing RHIC beam energy scan and planned FAIR experiments, one is motivated to reliably 
determine the thermodynamic functions at finite chemical potential. In this article we compute the NLO pressure of 
quarks and gluons at finite T and /i. The computation utilizes a high temperature expansion through fourth order in 
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the ratio of the chemical potential to temperature. This allows us to reliably access the region of high temperature and 
small chemical potential. We compare our final result for the NLO HTLpt pressure at finite temperature and chemical 
potential with state-of-the-art perturbative quantum chromodynamics (QCD) calculations and available lattice QCD 
results. 

The paper is organised as follows. In Sec. [Ill we will briefly review HTLpt. In Sec. Illll we discuss various quantities 
required to be calculated at finite chemical potential based on prior calculations of the NLO thermodynamic at zero 
chemical potential jsof . In Sec. II VI we reduce the sum of various diagrams to scalar sum-integrals. A high temperature 
expansion is made in Sec. |V]to obtain analytic expressions for both the LO and NLO thermodynamic potential. We 
then use this to compute the pressure in Sec I VII We conclude in Sec. IVIII Finally, in Appendices Rl and [B] we collect 
the various integrals and sum- integrals necessary to obtain the results presented in the main body of the text. 



II. HARD THERMAL LOOP PERTURBATION THEORY 



HTL perturbation theory [291 - 13 If is a reorganization of the perturbation series for hot and dense QCD which has 
the following Lagrangian density 

= {Cqcd + Chtl) ^ + A£htl , (8) 

where A£htl collects all necessary renormalization counterterms and £htl is the HTL effective Lagrangian [l^, . 
It can be written compactly as 

^HTL = -^(1 - S)mlTv (g^^ (jf^) ^'P^ + (1 - ^) *"^g^7" (^) ^ ' (9) 

where D is a covariant derivative operator, y = (l,y) is a light like vector and (• • •) is the average over all possible 
directions, y, of the loopmomenta. The HTL effective action is gauge invariant, nonlocal, and can generate all of the 
HTL n-point functions |22l. [23| . which are interrelated through Ward identities. The mass parameters rriD and niq are 
the Debye screening and quark masses in a hot and dense medium, respectively, which depend on the strong coupling 
g, temperature T, and the chemical potential ^. In the end we make expansions in mo and rtiq in ^ treating the 
masses as order g in order to make the calculation tractable (29l - [3ll | . The n"^ loop order in the HTLpt loop expansion 
is obtained by expanding the partition function through order and then taking (5 — > 1 [2^ - l35| . In this work, we 
will fix the parameters mo and mq by employing a variational prescription which requires that the first derivative of 
the thermodynamic potential with respect to both mo and mq vanishes such that the free energy is minimized. In the 
following, we generalize the NLO thermodynamic potential calculation from the case of zero chemical potential (soj 
to finite chemical potential. 



HI. INGREDIENTS FOR THE NLO THERMODYNAMIC POTENTIAL IN HTLPT 



The LO HTLpt thermodynamic potential, r^T o, for an SU{Nc) gauge theory with Nf massless quarks in the 
fundamental representation can be written as [2^, [SO] 

^LO = dAJ'g + dpTq + Ao^o , (10) 

where dp = NfN^ and rf^ = — 1 with Nc is the number colour. J-'q and J-'g are the one loop contributions to quar k 
and gluon free energies, respectively. The LO countcrterm is the same as in the case of zero chemical potential (29[ 

At NLO one must consider the diagrams shown in Fig. [51 The resulting NLO HTLpt thermodynamic potential can 
be written in the following general form (soj 

^^NLO = ^LO + dA [J-3g + J'4g + J'gh + J'gct] + dASp [^3qg + ^iqg] 

d d 

+dFJ'qct + Ai£:o + Aiml-^ni^o + Aim^ — f^LO , (12) 

where sp = Np/2. At NLO the terms that depend on the chemical potential are J-q, J-sqg, J'4qg, J^qt, Aim^, and 
Aim^ as displayed in Fig. [51 The other terms, e.g. J-'g, J-^g, J^ig-, J-gh and J-gct coming from gluon and ghost loops 
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FIG. 2. Diagrams containing fermionic lines relevant for NLO thermodynamics potential in HTLpt with finite chemical 
potential. Shaded circles indicate HTL n-point functions. 

remain the same as the /i = case (sOj- We also add that the vacuum energy counterterm, Aifo, remains the same as 
the /X = case whereas the mass counterterms, Aim|, and Aim^, have to be computed for /i 7^ 0. These counterterms 
are of order b. This completes a general description of contributions one needs to compute in order to determine NLO 
HTLpt thermodynamic potential at finite chemical potential. We now proceed to the scalarization of the necessary 
diagrams. 

IV. SCALARIZATION OF THE FERMIONIC DIAGRAMS 

The one-loop quark contribution coming from the first diagram in Fig. [5] can be written as 
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and Tp is defined by the following integral [30j 
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In three dimensions e — ?• and (fT6l) reduces to 



^iP^, iPp + IpI 
2|p|°^Po-|p| 



(17) 



with P — (Po,p). In practice, one must use the general form and only take the limit e — > after renormalization. 
The HTL quark counterterm at one-loop order can be rewritten from the second diagram in Fig. [2] as 
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The two-loop contributions coming from the third and fourth diagrams in Fig. [5] are given, respectively, by 

^3,g = ^ Tr [F^ (F, Q,R)S{Q)x (F, Q, R)S{R)] A^, (P) , (19) 



{PQ} 
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^^19 = ^.9' ^ Tr [r'-(P, -P, Q, g)5(g)] A^,(P) , (20) 

where S is the quark propagator and IS.^" is the gluon propagator. and F'^'' are HTL-resumnied 3- and 4-point 
functions. Details concerning the HTL n-point functions can be found in Refs. (29l. Isoj. 
In general covariant gauge, the sum of ((T9)) and ((20)) reduces to 

-^3,,+4,, - \9' ^ I Ax(P)Tr [r«°5(g)] - AT(F)Tr S{Q)V^ S{B!)\ 



{PQY 



-A^(p)Tr[r°5(g)r05(i?')] \ , (21) 



where A^ is the transverse gluon propagator, Ax is a combination of the longitudinal and transverse gluon propaga- 
tors [lOl, and R' = Q — P. After performing the traces of the 7-matrices one obtains [s^ 
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where Aq and are defined in ((T4|) and ([15]), respectively. We add that the exact evaluation of two-loop free energy 
could be performed numerically and would involve 5-dimensional integrations; however, one would need to be able 
to identify all divergences and regulate the numerical integration appropriately. Short of this, one can calculate the 
sum-integrals by expanding in a power series in m^i/T, m^/T, and [ijT in order to obtain semi-analytic expressions. 



V. HIGH TEMPERATURE EXPANSION 



As discussed above, we make an expansion of two-loop free energies in power series of rriD/T and viq/T to obtain 
a series which is nominally accurate to order . The HTL n-point functions can have both hard and soft momenta 
scales on each leg. At one-loop order the contributions can be classified "hard" or "soft" depending on whether the 
loop momenta are order T or gT, respectively; however, since the lowest fermonic Matsubara mode corresponds to 
Pq = ttT, fermion loops are always hard. The two-loop contributions to the thermodynamic potential can be grouped 
into hard- hard (hli), hard-soft (hs), and soft-soft (ss) contributions. However, we note that one of the momenta 
contributing is always hard since it corresponds to a fermionic loop and therefore there will be no two-loop soft-soft 
contribution. Below we calculate the various contributions to the sum-integrals presented in Sec. IIVI 



A. One-loop sum-integrals 

The one-loop sum-integrals (|T3|) and (|18p correspond to the first two diagrams in Fig. [2] They represent the leading- 
order quark contribution and order-(5 HTL counterterm. We will expand the sum-integrals through order taking 
ruq to be of (leading) order g. This gives a result which is nominally accurate (at one-loop) through order 5^.0 



^ Of course, this won't reproduce the full pQCD result in the limit g — )■ 0. In order to reproduce all known coefficients include , one 
would need perform a NNLO HTLpt calculation. 
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1. Hard Contribution 

The hard contribution to the one-loop quark self-energy in (jl3p can be expanded in powers of as 
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Note that the function Tp does not appear in term. The expressions for the sum-integrals in 
Appendix]^ Using those expressions, the hard contribution to the quark free energy becomes 
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Expanding the HTL quark counterterm in (|18p one can write 
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are listed in 
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where the expressions for various sum-integrals in ([25]) are listed in Appendix El Using those expressions, the hard 
contribution to the HTL quark counterterm becomes 



■^qct ~ 



(26) 



We note that the first term in cancels the order-e" term in the coefficient of in (P^ . There are no soft 

contributions either from the leading-order quark term in or from the HTL quark counterterm in (|18p . 



B. Two-loop sum-integrals 



Since the two-loop sum-integrals given in (j21[) contain an explicit factor of g^, we only require an expansion to 
order m'^mu/T^ and m^/T^ in order to determine all terms contributing through order g^. We note that the soft 
scales are given by rriq and mo whereas the hard scale is given by T, which leads to two different phase-space regions 
as discussed in Sec. IV Al In the hard-hard region, all three momenta P, Q, and R are hard whereas in the hard-soft 
region, two of the three momenta are hard and the other one is soft. 



The hh contribution 



The self-energies for hard momenta are suppressed [22, [23, [30] by nijj/T^ or rn^/T'^ relative to the propagators. 
For hard momenta, one just needs to expand in powers of gluon self-energies H^, Hl, and quark self-energy S. So, 
the hard-hard contribution of -Fsqg and Tiqg in (|2ip can be written as 
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+ 2my[d-l) 
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where the various sum-integrals are evaluated in Appendices [K\ and [Bl Using those sum-integral expressions, the hh 
contribution becomes 
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^. T/ie /is contribution 



Following Ref. |30j one can extract the hard-soft contribution from (|2ip as the momentum P is soft whereas momenta 
Q and i? are always hard. The function associated with the soft propagator Ar(0,p) or Ax(0,p) can be expanded 
in powers of the soft momentum p. For At(0, p), the resulting integrals over p are not associated with any scale and 
they vanish in dimensional regularization. The integration measure scales like the soft propagator Aj!f(0,p) 

scales like 1 / , and every power of p in the numerator scales like m £> . 

The contributions that survive only through order g^m^T and m^mug^T from J-3qg and -F^qg in (|2ip are 
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Using the sum- integrals contained in Appendices |X] and [B] the hard-soft contribution becomes 
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5. r/ie ss contribution 



As discussed earlier in Sec. |V] there is no soft-soft contribution from the diagrams in Fig. [2] involving since at least 
one of the loops is fermionic. 



C. Thermodynamic potential 



Now we can obtain the HTLpt thermodynamic potential r2(T, /x, a^, mu, m^, (5) through two-loop order for which 
the contributions involving quark lines are computed here whereas the ghost and gluon contributions are computed 
in Ref. (soj . We also follow the same prescription as in Ref. [soj to determine the mass parameter m/^ and rriq from 
respective gap equations but with finite quark chemical potential, /i. 
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1. Leading order thermodynamic potential 



Using the expressions of !Fq with finite quark chemical potential in (j24p and !Fg from Ref . [30| , the total contributions 
from the one-loop diagrams including all terms through order becomes 
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where m^i, m^, A. and fi are dimensionless variables: 
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Adding the counterterm in (jlip . we obtain the thermodynamic potential at leading order in the (5-expansion: 
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where we have kept terms of 0{t) since they will be needed for the two-loop renormalization. 



2. Next-to-leading order thermodynamic potential 



The complete expression for the next-to-leading order correction to the thermodynamic potential is the sum of 
the contributions from all two-loop diagrams, the quark and gluon counterterms, and renormalization counterterms. 
Adding the contributions of the two- loop diagrams, J-^qg+nqg^ involving a quark line in (|28p and pOI) and the contri- 
butions of T'ig+ig+gh boTCL Ref. [s^l , oue obtains 
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where ca = -/Vc and sp = Nf/2. 

The HTL gluon counterterm is the same as obtained at zero chemical potential [SC 
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The HTL quark counterterm as given by (pS)) is 
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The ultraviolet divergences that remain after adding ((37|) . (|38p . and ((39)l can be removed by renormalization of the 
vacuum energy density £q and the HTL mass parameter ttid and mq. The renormalization contributions (soj at first 
order in S are 
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Using the above counterterms, the complete contribution from the counterterms in (|40p at first order in 6 at finite 
chemical potential becomes 
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Adding the contributions from the two- loop diagrams in (|37p . the HTL gluon and quark counterterms in (|38p and (|39p . 
the contribution from vacuum and mass rcnormalizations in (j44|) . and the leading-order thermodynamic potential 
in p6p we obtain the complete expression for the QCD thermodynamic potential at next-to-leading order in HTLpt: 
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VI. PRESSURE 



In the previous section we have computed both LO and NLO thermodynamic potential in presence of quark chemical 
potential and temperature. All other thermodynamic quantities can be calculated using standard thermodynamic 
relations. The pressure is defined as 
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where rriD and niq are determined by requiring 
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This leads to the following two gap equations which will be solved numerically 
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For convenience and comparison with lattice data [6[, we define the pressure difference 

^P[T,^i) = P{T,y)-P{T,Q). 
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In Figs. [3] and m we present a comparison of NLO HTLpt pressure with that of four-loop pQCD [T^l as a function of 
the temperature for two and three loop running of as- The NLO HTLpt result differs from the pQCD result through 
order Inttg at low temperatures. A NNLO HTLpt calculation at finite fji would agree better with pQCD In as 
as found in /i = case [3ll |. The HTLpt result clearly indicates a modest improvement over pQCD in respect of 
convergence and sensitivity of the renormalisation scale. In Fig. [S] the pressure difference, AP, is also compared with 
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FIG. 3. The NLO HTLpt pressure scaled with ideal gas pressure plotted along with four-loop pQCD pressure [l3| for two 
different values of chemical potential with Nf — 3 and 2-loop running coupling constant as. The bands are obtained by varying 
the renormalisation scale by a factor of 2 around its central value A = ^-k^T"^ -f ji^ /-k^ We use A-j^ = 290 MeV based 

on recent lattice calculations [S^l of the three- loop running of as- 
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FIG. 4. Same as Fig. [3] but for 3- loop 



the same quantity computed using pQCD [l3| and lattice QCD @. Both LO and NLO HTLpt results are less sensitive 
to the choice of the renormalisation scale than the weak coupling results with the inclusion of successive orders of 
approximation. Comparison with available lattice QCD data Q suggests that HTLpt and pQCD cannot accurately 
account for the lattice QCD results below approximately STc; however, the results arc in very good qualitative 
agreement with the lattice QCD results without any fine tuning. 
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FIG. 5. (Left panel) AP for Nf = 3 is plotted as a function of T for two-loop HTLpt result along with those of four-loop 
pQCD up to In as and lattice QCD up to O (/i^) using 2-loop running coupling constant Us. (Right panel) Same as 
left panel but using 3-loop running coupling. In both cases three different values of /i are shown as specified in the legend. The 
bands in both HTLpt and pQCD are obtained by varying the renormalisation scale by a factor of 2 around its central value 
A = 27rVT2 + ^2/7r2 



VII. CONCLUSIONS AND OUTLOOK 

In this paper we have generalized the zero chemical potential NLO HTLpt calculation of the QCD thermodynamic 
potential |30| to finite chemical potential. We have obtained (semi-) analytic expressions for the thermodynamic 
potential at both LO and NLO in HTLpt. The results obtained are trustworthy at high temperatures and small 
chemical potential since we performed an expansion in the ratio of the chemical potential over the temperature. 

This calculation will be useful for the study of finite temperature and chemical potential QCD matter. This is 
important in view of the ongoing RHIC beam energy scan and proposed heavy-ion experiments at FAIR. Using the 
NLO HTLpt thermodynamic potential, we have obtained a variational solution for both mass parameters, uiq and 
rriD, and we have used this to obtain the pressure at finite temperature and chemical potential. When compared 
with the weak coupling expansion of QCD, the HTLpt pressure helps somewhat with the problem of oscillation of 
successive approximations found in pQCD. Furthermore, the scale variation of the NLO HTLpt result for pressure 
is smaller than that obtained with the weak coupling result. The HTLpt pressure shows some deviations from the 
lattice data below 3 Tc which suggests that the calculation should be extended to NNLO. In addition, getting better 
agreement with pQCD at low temperature will require going to NNLO. This is indeed a very challenging job which 
represents work in progress. 

We also note that, based on the results obtained herein, one can straightforwardly compute quark susceptibilities. 
In a forthcoming paper we will compare the NLO HTLpt results for quark susceptibilities with lattice data and other 
theoretical models of QCD matter. 
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Appendix A: Sum-Integrals 

In the imaginary-time (Euclidean time) formalism for the field theory of a hot and dense medium, the 4-momentum 
P = {Pq, p) is Euclidean with P^ = Pq +p^. The Euclidean energy Pq has discrete values: Pq = 2mrT for bosons and 
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FIG. 6. The contour corresponding to (|A.3P in complex Po plane. The crosses are poles of the thermal weight factor which are 
shifted by an amount /i from the Re[Po] axis. 



Po = i'2n + 1)ttT — for fermions, where n is an integer running from — oo to oo, fj. is the quark chemical potential, 
and T = 1//3 is the temperature of the medium. Loop diagrams usually then involve sums over Pq and integrals over 
p. In dimensional regularization, the integral over spatial momentum is generalized to d = 3 — 2e spatial dimensions. 
We define the dimensionally regularized sum-integral as 
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where 3 — 2e is the dimension of space, A is an arbitrary momentum scale, P is the bosonic loop momentum, and {P} 
is the fermionic loop momentum. The factor (e"' /Att)'^ is introduced so that, after minimal subtraction of the poles in 
e due to ultraviolet divergences, A coincides with the renormalization scale of the MS renormalization scheme. 

We describe below the technique of contour integration [l^ [l3| in the complex plane to evaluate the frequency sum 
over Pq. Consider a meromorphic function /(Po) that originates from a loop diagram, then one can write 



^ E /(P.)^T£g/(P„)f.a.,h 



Po={2n+l)7rT~ifj, 



^ \ ^' = X I X (2^*) E Residues , 



(A.3) 



provided f{Po) is regular in Re(iPo) ~ fi line as shown in Fig. [6] Below we demonstrate two examples, a simpler 
one involving only loop momentum and a complicated one involving fourth power of loop momentum and the HTL 
angular function, which would be relevant for evaluating sum-integrals: 

(i) Simpler one: 
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(ii) Involving HTL term: 
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After performing the frequency sum, one is left with dimensionaUy regularised spatial momentum integration, which 
are also discussed in Appendix |B] However, all other frequency sums can be evaluated in similar way as discussed 
above. 



1. Simple one loop sum-integrals 

The specific fermionic one-loop sum-integrals needed are 
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2. HTL one loop sum-integrals 



We also need some more difficult one-loop sum-integrals that involve the HTL function defined in 
The specific fermionic sum-integrals needed are 
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3. Simple two loop sum-integrals 
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The generalized two loop sum-integrals can be written from [30j as 
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After applying Eq. (jA.27[) and using the delta function to calculate the Pq and Qo integrations, the sum-integral 
(|A.23P reduces to 
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and using the result of Eq. (|B.10[) . we get sum-integral (jA.19[) and agree with [l7| . 
After applying Eq. (|A.27p . the sum- integral (|A.20p reduces to 
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where np{p) = np{p) + n^{p) . Now using the result of 4-dimensional integrals from [s^l and applying Eq. (|B.3[) and 
Eq. (|B.5p . we can calculate sum- integral Eq. (jA.20|) . The sum- integrals (|A.21[) can be calculated in same way: 
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The sum-integral (jA.22p can be written as 
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Using Eq. (|A.27p and after doing Pq and Qq integrations, first sum-integral above reduces to 



PqQo 

p2Q2^4 = 

{PQ} pq 



?ijp) - npip) np{q)-n+{q) pq 



2p 



2q 



(A.31) 



(A.32) 



(A.33) 



and the result is given in Eq. (|B.9I) . The second term and third terms sum- integrals above are linear combinations of 
Eq. (jAT20l) and Eq. (|X2T|) . Ad ding a ll of them, we get required sum-integral. 
Similarly after applying Eq. (jA.27p . the sum- integral (|A.23p reduces to 
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Using the 4-dimensional integrals from [30j | and Eqs. (jB.2p . (|B.3p . (|B.4p and (jB.6p . we obtain the sum- integral (|X23t . 
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Similarly after applying Eq. (jA.27|) . the sum- integral (|A.24[) reduees to 
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Using the above angular average, we find 
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Similarly after applying Eq. (|A.27p . the sum- integral (jA.25[) reduces to 



{PQ} 



p^P^Q^K^ 



nsip) 



Q^r^R^ 



"-f(p) 
2p 



Q^K^ \p^ g2 



Po = ~ip 



nF{p)nFiq) "T^ — p^ — q^ 



Apq 



p2 A{p,q,r) 
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(A.42) 



Now 
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(A.43) 



and 



2 9 
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9 9 9 

— p'^ — q^ 
A{p,q,r) 



1/^ 

2e \r4 



(A.44) 



Using the above angular average, we have 
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(Ib31). 



nF{p)nB{q) if_ 
2pq r4 



Po =—ip 



(A.45) 



Using the 4-dimensional integrals from [30] and Eqs. (jB.2[) . (|B.3[) . (|B.4p and (jB.Sp . we obtain the sum- integral (jA.24l) 
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4. HTL two loop sum-integrals 
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(A.46) 
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Appendix B: Integrals 

1. Three dimensional integrals 

We require one integral that does not involve the Bose-Einstein distribution function. The momentum scale in 
these integrals is set by the mass m = mjj. The one- loop integral is 
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2. Thermal Integrals 
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Thermal integrals containing the triangle function: 
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Thermal integrals containing both the triangle function and HTL average are listed below: 
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Second set of integrals involve the variables J'c = |p + ^./c\'■ 
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The integral (|B.12p can be evaluated directly in three dimensions at finite chemical potential. The other integrals 
Eqs. ()B.13p - (|B.20P can be evaluated following the same procedure as discussed in [l^l at finite chemical potential. 
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